7-ray burst after-glow: Confirming the cosmological fireball model 
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ABSTRACT 

The recent detection of delayed X-ray and optical emission, "afterglow," associated 
with 7-ray bursts (GRBs) supports models, where the bursts are produced by 
relativistic expanding blastwaves, "fireballs," at cosmological distances. The detection 
of absorption lines in the optical afterglow of the GRB of 8 May 1997 confirms that 
the sources lie at cosmological distance. We show here that the new features detected 
in GRB970508 afterglow, radio emission one week following the burst and a 2 day 
increase in optical flux, are consistent with the blastwave model. The fireball optical 
depth at radio frequencies is much smaller than previously estimated, which accounts 
for the observed radio emission. The initial suppression of optical flux is consistent 
with that predicted due to electron cooling. The combined radio and optical data 
imply that the fireball energy is ~ lO^'^erg, and that the density of the medium into 
which the blastwave expands is ~ lcm~'^, a value typical for gas within galaxies. We 
predict the time dependence of the radio flux and the absorption frequency, which 
constitute tests of the fireball model as described in this paper. 
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1. Introduction 

The origin of GRBs, bursts of 0.1 — IMeV photons lasting for a few seconds, remained 
unknown for over 20 years ( Fishman Sz Meegan 199"5| ), primarily because GRBs were not detected 
until this year at wavebands other than 7-rays. Phenomenological considerations based on 7-ray 
data were used to argue that the bursts are produced by relativistic expanding blastwaves at 
cosmological distances ( Paczyiiski 1986| , Goodman 1986| , Rees &i Meszaros 1992, Piran 1996). The 



fireball model, which predicts delayed emission at wavelengths longer than 7-rays (Paczyhski 
fc Rhoads 1993| , [Katz 199l| , [Meszaros Sz Rees 1997] , [Vietri 1997| ), has gained support ( [Waxmari 
1997, Wijers, Rees fc Meszaros 1997| ) from the recent detection (van Paradijs et al. 1997) of delayed 



X-ray and optical emission associated with GRBs. The key new development is the availability 
from the BeppoSAX satellite of accurate positions for GRBs shortly after their detection. 



GRB970508 was detected by the BeppoSAX satellite ( |Piro, Scarsi fc Butler l"995| ) on 8 May 
1997. The burst lasted for ~ 15s, with 7-ray fluence ~ 3 x 10~^ergcm~^ carried mainly by photons 
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of energy ~ 0.5MeV ( [Kouveliotou at al. 1997 ). Following the detection in 7-rays, X-ray ( piro et 
al. 1997| ), optical (Bond 1997, Djorgovski et al. 1997a, pjorgovski et al. 1997b , Djorgovski 
al. 1997c| , [Mignoli et al. 199% phevalier fc Ilovaisky T997| ), and radio ( [Frail &: Kulkarni 199% ) 



emission varying on time scale of days was observed from the direction of the GRB. The location 
of the GRB is determined with accuracy better than ~ 3', while the locations of the X-ray, optical 
and radio sources are determined to better than 50", 1", and Imas respectively. The fact that the 
locations of all sources coincide, and the unusual variability observed following the GRB, strongly 
suggest that the X-ray, optical and radio emission are associated with the object producing the 
GRB. A spectrum of the optical transient, taken 2 days after the GRB, shows a set of absorption 
features, associated with Fe II and Mgll and shifted to long wavelength, implying that the 



absorbing system lies at a cosmological redshift z=0.835 (Metzger et al. 1997). This sets a lower 
limit to the GRB redshift, z > 0.835, and to the energy emitted by the source in 7-rays (assuming 
isotropic emission), > 10^^ erg. Intensive monitoring of GRB970508 revealed two new afterglow 
features: radio emission was observed one week following the burst ( Frail &: Kulkarni 1997| ), 



and the optical flux was observed to increase for 2 days following the burst ( Djorgovski et al] 
1997a| , pjorgovski et al. 1997b| ). 



The early detection of radio emission, combined with optical data, appears to be inconsistent 
with published fireball model predictions ( Paczyiiski fc Rhoads 199"3| , [Meszaros Sz Rees 1997 ). 
We show here, however, that the fireball optical depth at radio frequency is significantly smaller 
than previously estimated, and that the same process responsible for the optical afterglow may 
also produce the observed radio emission. The combined radio and optical data are shown to be 
consistent with the fireball model, and to provide information on the fireball parameters and on 
the ambient medium into which the blastwave expands. The suppression of optical flux at early 
times provides evidence for the predicted suppression due to electron cooling ( Waxman 1997| ), and 
suggests that the GeV emission observed for several hours following several strong GRBs ( [Hurley 
et al. 1994) is produced by inverse-Compton scattering of afterglow X-ray photons. 



Our goal is to present a qualitative analysis of the new afterglow features observed in 
GRB970508, in order to identify the physical processes that are responsible for the observed 
behavior. We therefore adopt a simple approximate description of the fireball expansion, that 
permits the derivation of the main results without complicated calculations. Our values for 
numerical parameters should be considered as order of magnitude estimates. 



2. The fireball model and predicted X-ray, optical and radio afterglow 

The underlying source, which produces the initial explosion that drives the expanding fireball 
is unknown (although several plausible candidates have been proposed, see, e.g., Piran 1996 for 
review). However, the 7-ray observations suggest the following scenario for the emission of the 
observed 7-rays ( Paczyhski 1986| , [Goodman 1986 , Rees Sz Meszaros 1992). A compact, tq ~ 10''cm, 
source releases an energy E comparable to that observed in 7-rays, E ~ 10^^ erg, over a time 
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T < 100s. The large energy density in the source results in an optically thick plasma that expands 
and accelerates to relativistic velocity. After an initial acceleration phase, the fireball energy is 
converted to proton kinetic energy. A cold shell of thickness cT is formed and continues to expand 
with time independent Lorentz factor 7 ~ 300. The GRB is produced once the kinetic energy 



is dissipated at large radius, r > 10 cm, due to internal collisions within the ejecta (Paczyhski 



fc Xu 1994 , Meszaros &: Rees 1994 ) and radiated as 7-rays through synchrotron and possibly 



inverse-Compton emission of shock accelerated electrons. 

Following internal collisions, which convert part of the energy to radiation and which result 
from variations in 7 across the expanding shell, the fireball rapidly cools and continues to expand 
with approximately uniform Lorentz factor 7. As the cold shell expands it drives a relativistic 
shock (blastwave) into the surrounding gas, e.g. into the inter-stellar medium (ISM) gas if the 
explosion occurs within a galaxy. In what follows, we refer to the surrounding gas as "ISM gas," 
although the gas need not necessarily be inter-stellar. The shock propagates with Lorentz factor 
7s = 2^/^7, and behind it the (rest frame) number and energy densities of the shock heated ISM 
are n' = 4771 and e' = 47^nmpC^, respectively, where n is the ISM number density ahead of the 
shock. The width of the shock heated ISM shell is r/47^, where r is the fireball radius. 

At late times, the blast wave approaches a self similar behavior. The expansion is well 
approximated by the self-similar solution for radii r > Tc, where ( Waxman 1997| ) 



^ 2,<lQ^^(^^\'\m. (1) 
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Here the fireball energy is E = W^^E52erg, the ISM density is n = nicm ^, and T = IOTiqs. We 
have chosen to give our numeric results using a fireball energy value that is somewhat higher than 
the typical energy observed in 7-rays, ~ lO^^erg, since the conversion of fireball energy to 
7-rays is not expected to be 100% efficient. For r > rc, the shell Lorentz factor is 



7 



^ ^ (r/r.)-'/' = 190 f ) ' (r/r,)-3/2. (2) 
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Photons emitted from the shell at radius r are seen by a distant observer at a time 
t = r/2j'^c = T{r/rc)'^ after the GRB, with arrival time spread comparable to t. The main 
contribution to the time delay and spread is due to two effects. First, the radiation seen by a 
distant observer is emitted from a cone of the fireball around the source-observer line of sight, with 
an opening angle ^ I/7. Photons emitted from such a cone are spread over t = r/2^'^c. Second, 
as we show below, during most of the afterglow the synchrotron cooling time of electrons is larger 
than the fireball expansion time, ~ r/jc as measured in the fireball frame. Thus, electrons emit 
radiation over a time ~ r/jc in the fireball frame, corresponding to a time ~ r/2^'^c as seen by a 
distant observer (Note, that the delay due to the difference ts between the shock propagation time 
to radius r and r/c, ts = r/167^c, is not the main factor determining the time t at which radiation 
from radius r is observed). 
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The shock driven into the ISM continuously heats new gas, and produces relativistic electrons 
that may produce the delayed radiation observed on time scales of days to months. In order to 
calculate the synchrotron emission from the heated ISM shell we need to determine the magnetic 
field and electron energy within the heated shell. We assume that the magnetic field energy density 
(in the shell rest frame) is a fraction of the equipartition value, -B^/Svr = ^se', and that the 
electrons carry a fraction of the energy. Since the Lorentz factor associated with the thermal 
motion of protons in the shell rest frame is 7, this implies that the Lorentz factor of the random 
motion of a typical electron in the shell rest frame is jem = iel'fnp/'^e- With these assumptions, 
and using eqs. (||) and (|2|), we find that the observed frequency of synchrotron emission from 
typical electrons, Um — 77em^-^V2^"^eC, is 

1 /2 

where t = Ifdaydays, and z is the cosmological redshift of the burst. We have chosen to give 
numerical results using ~ ■Cs ~ 0.1, since such values are typically required for the production 
of the GRB itself. We show below that the electron synchrotron cooling time is longer than the 
dynamical time, td = r/jc, the time for significant fireball expansion. In this case, the observed 
intensity at is (Waxman 1997) 
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(CB/0.1)i/2i?52mJy. (4) 



Here (and throughout the paper) we assume a flat universe with Hubble constant Hq = 75km/s Mpc. 
The Jy flux unit is IJy = 10~^^ erg/cm^ sHz. 

A natural prediction of the fireball model, see eqs. (^) and (Q), is optical emission at a level 
of ImJy at a delay of order Iday following the GRB. 

Emission at > is produced by electrons with Lorentz factor higher than 7em- 

If the 

electron distribution follows a power-law, dNe/dje oc 7"^ for 7e > ^em, as expected for shock 
acceleration, then for v > Vm 

F, = F,Jv/v^{t)\-^, (5) 

with a = {p — l)/2. Although p is expected to be similar for the GRB and for the afterglow, 
a during the after-glow is expected to be smaller by 1/2 due to increase in cooling time: the 
synchrotron cooling time must be short compared to the dynamical time during the GRB, resulting 
in a = p/2, while it is long during the afterglow, giving a = {p - l)/2 ( IWaxman 19"9^ ). For the 
GRB a ~ 1, implying a ~ 0.5 for the after-glow. 

The synchrotron emission of electrons with Lorentz factor 7em is concentrated mainly at 
frequencies v ^ Vrn- However, the emission extends to lower frequencies, with power radiated 
per unit frequency proportional to {v /vmY^^ ■ Thus, we expect the synchrotron flux to extend to 
u < Urn following eq. (|5|) with a = —1/3. This implies that the flux at a fixed frequency v increases 
with time as t^^^ as long as v < Vm- We have so far assumed that the fireball optical depth to 
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synchrotron absorption is small. This may not be the case for low frequencies. The synchrotron 
optical depth at is small, = E^-^^'^nit^J^. For > z/^, t oc z/~(p+^)/^, while 

T oc for 1/ < i/fn due to the u^^^ low frequency synchrotron tail. The frequency for which the 

3/5 

optical depth is 1, i^a = i^mTm , is 

= 1 ' iCe/0.2)-\^B /O.iy /^Eli'nl^' GHz, (6) 

and the flux at z^^ < ^ J^m, Fa = 2-Fjy„(z^/t'm)^^^, is 



= 0.3 



1-1/^2 

1- 



(ee/0.2)-2/3(^^/0.1)l/3^5/6^1/2^^^^^^^j^^l/2^^^^Q(^jj^^l/3 ^jy_ 



(7) 

Eq. (|6|) and (0) indicate, that radio emission at a level of ImJy is expected at the ~ lOGHz range 
on time scale of weeks. The emission should be suppressed below ~ IGHz due to high optical 
depth. 

We note here, that a population of electrons with 7e <C 7em may also contribute to the 
flux and optical depth at v <^ Vm- Although in the model described in this paper the electron 
distribution is dominated by electrons with Lorentz factor 7e ~ 7em, the distribution may extend 
to low energy, 7e <C Jem, without significantly affecting the results provided that low energy 

electrons constitute a small fraction of the electron population. The distribution may extend, e.g., 
to 7e <C 7em as dNe/d'je OC j^'''' with < p' ^ 1 or < (for p' ^ 1 the fraction of electrons with 
7e ~ lem is Small). While the contribution of such low energy electron population to the flux at 
^ < ^m, Fy oc is not large compared to that given in (|^, their contribution to the optical 

depth at < z^rra; oc may result in self-absorption frequency which is significantly 

higher than (^). For p' = we have va = VmTm'^ , i.e. 

i^A = 4 [^y'^' fe/0.2)-V2(^^/0.l)V4ii;V4„V2(i/i^eek)-V4 GHz. (8) 

Eq. (|8|) implies, that the existence of low energy electron population does not change the 
conclusion that radio emission at a level of ImJy is expected at the ~ lOGHz range on time scale 



of weeks. (Our estimate for ua is significantly lower than previous estimates [e.g. Paczyiiski & 



Rhoads 1993 , Meszaros &: Rees 1997 | since we have taken into account the fact that the electron 



distribution flattens below 7e 



3. GRB970508 

Let us now compare the fireball model predictions with observations of GRB970508. The 
detected absorption lines imply that the GRB source redshift is z > 0.835, and the absence of 



Lyman-a lines imply z < 2.1 ( Metzger et al. 1997 ). The absence of CIV absorption further implies 



- 6 - 



z < 1.6 (N. Arav and D. Hogg, private communication). We therefore adopt z = 1 as the GRB 
source redshift. 

The radio flux detected 6 days after the GRB is weh described by 0.2(z^/3GHz)mJy in the 
range of 1 to lOGHz (Frail & Kulkarni 1997). The i^^ dependence, which is steeper than z^^/^ 



expected from the low frequency tail of synchrotron emission in the absence of absorption, implies 
that the radio emission is partly absorbed, with stronger absorption at lower frequencies due to 
increase in optical depth. In our model, oc z/^/^ for v/v^ ^ 1, and Fy oc for iz/i/^i <C 1 (in the 
presence of low energy electrons, dNe/dje oc 7"^' with p' > —2/3 for 7e ^ Jem, the dependence 
of flux on frequency is stronger, F^, oc with x > 2 for u/ua <^ 1). Radio observations therefore 
indicate that IGHz < va < lOGHz. The observed self-absorption frequency and the level of 
detected radio flux are in agreement with model predictions (P-^)- The upper limit of 7mJy at 6 
day delay at 90GHz (Shepherd, Metzger & Kulkarni 1997) is also consistent with the model i/^/^ 



dependence at > ua, which predicts a flux at 90GHz of ImJy. At delays shorter than 6 days, 
only upper limits to the flux at 1.4GHz are available. The upper limits are comparable to the 
flux detected at 6day delay, and therefore, while being consistent with the model, do not allow 
to establish the transient nature of the radio emission or test the model prediction, that the flux 
should be increasing with time. 



The R band, = 4 x 10 Hz, observations at 2.0 to 5.5 day delay ( Djorgovski et al. 



1997a| , pjorgovski et al. 1997"b| , pjorgovski et al. 19971 , |Mignoh et al. 199% [Chevalier fc Ilovaisky 



1997) are well described by a power- law, Fr = 37;uJy(t/2day)-i-3. The power law behavior is 
consistent with the model prediction (|5|). The approximately constant flux during 1 to 2 days 
delay indicate that the peak of the synchrotron emission Vrn passed through the R band at a 
delay tR ~ Iday. This is consistent with the model prediction (^). (Note that extrapolation 
of the t~^''^ behavior to t = Iday results in a flux ~ 2 times higher than observed. However, 
since we do not expect a sharp peak with in the flux observed at a given frequency as Vm drops 
below this frequency, due to the spread in photon arrival times, the approximately constant flux 
observed near the peak on time scale comparable to the delay is consistent with the model). The 
normalization of the power-law fit to the R band observations is ~ 3 times lower than implied by 
6q- (@)i d), and (|5|). However, given the simple description of the fireball behavior adopted in this 
paper, one should not draw conclusions based on this numerical discrepancy. The t~^'^ decline of 
R band fiux implies, through (|5|), a ~ 0.8, or p = 2.6. The implied frequency dependence of the 
flux is consistent with observations at other optical bands, although the narrow frequency range 
does not allow accurate determination of a from the frequency dependence. 



Consider next the optical ( [Bond 1997| , Djorgovski et al. 1997a) detection at ~ 7hr delay, 



and the X-ray detection at lOhr delay ( Piro et al. 1997| ). At early time, t < Iday, the R band 



frequency is smaller than i^m, the frequency of radiation emitted by the typical fireball electrons, 
and the R band flux should increase as t^/^. Using the the power-law fit to the observed R-band 
flux for t > Iday, and a t^/^ scaling for t < Iday, the predicted flux at 7hr delay is 50/iJy, 
signiflcantly above the observed flux, 11.1 it 2.2/iJy. This discrepancy may indicate the detection 
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of the predicted (Waxman 1997) suppression of optical flux at early times due to rapid electron 
cooling. The ratio of synchrotron cooling time, tg = GTrrriec/aT^eB''^ , to fireball deceleration time, 
td ~ 7/77 = ^r/Sjc, is tg/ta ~ 2(^e/0-2)~^(.^B/0.1)~^ii^52^''^*d{,y- Thus, our assumption that the 
synchrotron cooling time is longer than the dynamical time is valid for t ^ Iday. At 7hr delay 
the cooling time is comparable to the deceleration time, and the electrons cool on a dynamical 
time scale. At this time, synchrotron emission is significantly suppressed due to inverse-Compton 
emission ( Waxman 1997]) , by which electrons lose a significant fraction of their energy to the 
production of photons of high energy, exceeding 1 GeV during the first hours following the burst. 
This may account for the delayed GeV emission observed in several strong bursts. Finally, the 
X-ray emission detected at ~ lOhr delay ( Piro et al. 1997] ), with a flux of ~ 0.03/iJy at lO^^Hz, 
is consistent with the above model when the suppression due to electron cooling is taken into 
account. 



4. Conclusions 

The fireball model for GRBs is consistent with afterglow observation in X-ray, optical and 
radio wave bands, provided the kinetic fireball energy is E ^ lO^^erg, the density into which the 
blastwave expands is ~ lcm~^, and provided that the fraction of energy carried by electrons and 
magnetic field is .^e ~ 0.1. The inferred kinetic energy is consistent with the observed 7-ray 
energy, ~ lO^^erg, and implies that the efficiency with which kinetic energy is converted to 
7-rays is of order 10%. The density of the ambient medium is typical for interstellar gas, therefore 
suggesting that the explosions take place within galaxies. The requirement for significant energy 
to be carried by electrons and magnetic field is consistent with what is usually assumed for the 
production of the GRB itself. We have implicitly assumed that the fireball is spherical. However, 
our analysis is valid also for the case that the fireball is a jet, as long as the jet opening angle 
is larger than I/7. The fireball Lorentz factor during the radio observations is not very large, 
7 = 4(£'52/ni)"'^/^(t/lweek)~^/^. This implies that the opening angle of the jet is not very small. 

Our model predicts that the frequency, at which the optical depth for synchrotron absorption 
is unity, should change slowly with time, i^a oc with < 2; < 1/4, and that on time scale of 
weeks IGHz < i^^ < lOGHz (cf. eq. The radio flux at v > ua should rise as t^^"^ following 

the GRB. These predictions can be tested with more frequent radio observations following future 
GRBs. (For GRB970508 the observed self-absorption frequency and radio flux at 6day delay are 
consistent with model predictions. However, at shorter delays only upper limits are available, 
which do not allow to establish the transient nature of the source or test the above predictions). 
Since the Lorentz factor after one week is not large, deviations from the simple scaling laws derived 
here may be observed at later times, due to the deceleration of the fireball to velocity which is not 
highly relativistic. 
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